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The rapid progress in precisely designing the surface dec- 
oration of patchy colloidal particles offers a new, yet un- 
experienced freedom to create building entities for larger, 
more complex structures in soft matter systems. However, 
it is extremely difficult to predict the large variety of or- 
dered equilibrium structures that these particles are able 
to undergo under the variation of external parameters, 
such as temperature or pressure. Here we show that, by 
a novel combination of two theoretical tools, it is indeed 
possible to predict the self-assembly scenario of patchy col- 
loidal particles: on one hand, a reliable and efficient op- 
timization tool based on ideas of evolutionary algorithms 
i helps to identify the ordered equilibrium structures to be 
expected at T = 0; on the other hand, suitable simulation 
techniques allow to estimate via free energy calculations 
the phase diagram at finite temperature. With these pow- 
erful approaches we are able to identify the broad variety 
of emerging self-assembly scenarios for spherical colloids 
decorated by four patches and we investigate and discuss 
the stability of the crystal structures on modifying in a con- 
trolled way the tetrahedral arrangement of the patches. 

J Experimental and theoretical investigations have provided un- 
ambiguous evidence that colloids with chemically or physi- 
cally patterned surfaces (commonly known as "patchy parti- 
cles") are very promising mesoscopic entities that can be used 
in hierarchical self-assembly processes to build up colloidal 
super- structures^. The anisotropy in the interactions of such 
particles in combination with the limited functionality and 
selectivity of the bonds offer unlimited possibilities for self- 
assembly scenarios. Thus, patchy particles are celebrated "to 
become the elementary brick of tomorrow's self-assembled 
materials" 3 , with promising applications in photonic crystals, 
drug-delivery, electronics 4 , biomaterials, or catalysis^. 

The possibility to tailor the interactions of patchy parti- 
cles almost deliberately represents the basis for "bottom-up" 
strategies which allow to build up materials with desired prop- 
erties, starting from adequately designed units. Suitable ex- 
perimental techniques^"- allow to position patches on the col- 
loidal surface and to define their spatial extent with high pre- 
cision. A very impressive example for such an approach is a 
recent work on triblock Janus particles 5 : after decorating col- 
loids with two hydrophobic caps of tunable area, particles self- 
organize in the two-dimensional Kagome lattice target struc- 
ture. Complementary computer simulations^*^ have provided 
a complete phase diagram of the system, including also the 
disordered, fluid phase. 

For the case of triblock Janus particles the self-assembly 
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scenarios were easy to "guess". However, for more complex 
decorations and three dimensional systems, it is considerably 
more difficult to identify all ordered structures into which the 
particles might self-assemble. Semi-empirical approaches, ap- 
plied over many years in hard matter physics, rely on a pre- 
selection of candidate structures, based on experience, intu- 
ition or plausible arguments. In view of the rich wealth of 
unexpected ordered structures in soft matter systems such a 
procedure is bound to fail. In this contribution we propose 
a novel approach which helps to predict with high reliability 
the ordered equilibrium structures of a system of patchy col- 
loidal particles. We evaluate the phase diagram for different 
tetrahedral patch arrangements based on a highly reliable pre- 
selection step, and we investigate the stability of the emerging 
structures. We provide a comprehensive classification of these 
structures in dependence of the patch geometry and interpret 
them in terms of the strong competition between energetic, 
entropic and packing contributions. 

We use a standard model for patchy particles, introduced by 
Doye et al.—, which has been used in numerous studies (e.g., 
Refs.— ~— ). An isotropic Lennard- Jones (LJ) potential (speci- 
fying the spherical colloids) is modulated by an orientationally 
dependent factor; the latter mimics the patches which can be 
"located" on arbitrary positions and with variable extent on the 
surface of the spherical colloidal particle. We use the LJ pa- 
rameters o and 8 as units for length and energy. As in previous 
work—, the patch width is fixed by the parameter G pw = 0.3, 
leading to rather narrow patches - the interaction strength of 
a bond decreases to half its value when the patch alignment 
deviates by an angle of 20 degrees from the perfect one. Ad- 
ditionally, the potential is truncated and smoothed at a cutoff 
distance of r cut = 2.5a in order to prevent long-range interac- 
tions that are very unrealistic for colloids. Usually colloidal 
interactions are even of much shorter range, but this choice al- 
lows us to perform a comparison with previous results, if avail- 
able. We also note, that the long-range contributions (which 
we define by r > 1.5 

^min ^ l-68o, where r m j n is the minimum 
of the LJ potential) to the binding energies of the structures 
discussed in the following are well below ten percent of the 
total binding energy and of the same order of magnitude in all 
cases. In this contribution we consider only one value of the 
interaction range; however it has been shown, that a change 
in the interaction range can affect the regions of stability of 
the emerging phase s 1 ^ 20 . Packing fractions T| are calculated 
as r| = (n/6)(Na 3 /V), V and N being the volume of and the 
number of particles in the primitive cell, respectively. Bind- 
ing energy and pressure values are calculated as dimensionless 
quantities (£/* = U/ (Ne) and P* = Pc 3 /e, respectively). 

Here we focus on four-patch particles: introducing a geo- 
metrical parameter g, we vary the patch positions on the sur- 
face, ranging from a rather elongated to a flat, compressed 
tetrahedral arrangement. More specifically, the first patch pi 



is fixed at the north pole of the particle, while the three re- 
maining patches (with k G {2, 3, 4}) are located at a degree 
of latitude related to g: 
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See Figure Q] for a visualization. 
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Fig. 1 (Colour online) Visualization of the decoration of the patchy 
particles (large spheres) by the four patches (small spheres) as the 
geometric parameter g varies. Here (and in the following figures) 
blue spheres specify the patch located on the north pole of each 
particle (pi), while green spheres specify the remaining three 
patches (p^). The arrow represents the g value for which the patches 
form a regular tetrahedron. 

Our investigations display for the first time the surprisingly 
broad variety of ordered equilibrium structures that the parti- 
cles are able to form when the regular tetrahedral arrangement 
of the patches is modified via the parameter g. 

In order to investigate the self-assembly scenarios of the 
system, we use a two-step approach, combining two efficient 
and reliable numerical tools. First, we identify at vanishing 
temperature T possible ordered equilibrium structures by ap- 
plying an optimization tool which is based on ideas of evo- 
lutionary algorithms (EAs) and searches essentially among all 
possible lattice structures 21 . Second, we use suitably designed 
computer simulation s 13 ^ 23 to evaluate the free energy of the 
candidate structures previously identified, as well as the free 
energy of disordered phases, leading to the phase diagram at 
finite T. By performing full free energy calculations, we are 
thus including entropic contributions at nonvanishing temper- 
ature. 

In numerous applications to a wide variety of systems, EAs 
have turned out to be reliable, efficient and robust^2i~— : they 
cope well with high dimensional search spaces and cost func- 
tions that depend on the search space coordinates in highly 
nontrivial ways. We note that a different structure prediction 
technique, based on Monte Carlo (MC) simulations, has been 
introduced in the recent literature^. For a thorough compari- 
son between these two approaches, applied to patchy systems, 
we refer to Ref.—. 

Working at constant pressure P, the lattice parameters, the 
positions of the (up to eight 0) particles within the primitive 



cell and their orientations are optimized by minimizing the 
Gibbs free energy, which reduces to the enthalpy at T = 0. To 
cope with the large number of parameters, we use a phenotype 
implementation of EA based optimization technique s 1 ^ 25 . In 
short, such a scheme combines evolutionary global optimiza- 
tion with local optimization steps. In the former, a set of can- 
didate solutions (a "population" consisting of "individuals"), 
each characterized by its search space coordinates and a fitness 
value (the latter being a function of its Gibbs free energy), is 
undergoing an evolutionary process; recombination and mu- 
tation operations are used in order to create new individuals, 
whose fitness values are expected to gradually improve. The 
local optimization is carried out by a quasi-Newton algorithm, 
namely L-BFGS-B 29 -^. The algorithm described in Ref. 18 
has been augmented in the following way: (i) we make use of 
the angle-axis description for handling rotations of rigid bod- 
ies first introduced for soft matter systems with anisotropic 
particles in Ref. 31 . A detailed description and examples for 
successful application can be found in Ref s. 2221. ({{) Since 
the time-consuming local optimization steps are independent 
from each other, they can be carried out simultaneously on 
different processors. Further, we drop the concept of genera- 
tions in the algorithm and use an implementation of the pool 
based approach described in Ref.— instead, (iii) Along the 
evolutionary process, bond order parameter s 34 ^ 5 are used to 
distinguish between energetically equivalent, but structurally 
different lattices in order to retain a diverse population by em- 
ploying structural niches 36 . The algorithm does not only keep 
track of the global minimum structure, but also of low-lying 
local minima, which might become dominant at finite T— . 
As in any nondeterministic optimization approach there is no 
guarantee that the lowest identified minimum is indeed the 
global one, therefore, in a strict sense, we have to consider 
our structures putative groundstates^S. However, we affirmed 
our results by performing a large number of independent opti- 
mization cycles at each state point. 

The structures suggested by the EA are replicated as starting 
configurations for NPT MC simulations, with typically about 
500 particles in the simulation box. The equations of state of 
the fluid and solid phases are calculated using these simula- 
tions. For solids with non cubic symmetry the shape of the 
simulation cell is allowed to change during the simulations to 
obtain the equilibrium structure of the solid at each thermody- 
namic state. Free energies of solids are calculated using the 
Einstein molecule approac h 22 i 23 i 39 , based on modifications of 
the Frenkel-Ladd method 40 . The orientational degrees of free- 
dom need special treatment, explained in further detail in the 
supplementary material^. The free energy of the liquid is ob- 
tained by thermodynamic integration from the ideal gas. Once 
the free energy at a state point is known, free energies at other 
states are calculated by thermodynamic integration. Coexis- 
tence points are obtained by finding the temperature and pres- 
sure at which the chemical potential of the two phases in co- 
existence are equal. Once a coexistence point is known, the 
whole coexistence line can be integrated by using the Gibbs- 
Duhem method—. 

We start the discussion of our results with the ordered equi- 



^For higher particle numbers, the computational cost of the optimizations 



rapidly increases. All the structures presented in this contribution can be real- 
ized with no more than four basis particles. 




Fig. 2 (Colour online) Contour plot of the reduced energy U* (left) and packing fraction r| (right) as functions of the dimensionless pressure 
P* and the geometric factor g; colour codes are displayed (for r|, the range is cut off at r| = 0.45, i.e., all values below 0.45 appear in the same 
colour). White boundaries indicate the limits of stability of the respective ordered structures on the underlying grid, "be" stands for 

body-centered lattices, "fc" for face-centered lattices. Solid arrows specify the g value corresponding to a regular tetrahedral patch 
arrangement, dotted arrows correspond to g = 127.50, for which we present results at finite T. Labels specify configurations referred to in the 
text. 



librium structures at T = C0. The broad variety of identified 
lattice structures is summarized in Figure [2 The topologies 
of the energy- and packing fraction-landscapes are the result 
of a complex competition between two mechanisms: packing 
(minimizing V) vs. bond saturation (minimizing the energy 
U*). Accordingly, we can identify three regions in (P*,g) 
space: 

(i) For pressure values up to P* ~ 4.00 fully bonded struc- 
tures (i.e., U* ~ —2.00) dominate, which have rather small 
packing fractions, from r\ ~ 0.27 for open structures to r\ ~ 
0.64 for distorted body-centered (be) structures (compared to 
r| c± 0.74 for close-packed spheres). The spectrum of identi- 
fied lattices ranges from open, layered structures (g ~ 90.00) 
over be lattices (which dominate over a broad g range, i.e., 
93.75 < g < 135.00) to layers of hexagonally arranged parti- 
cles for g > 135.00. (ii) For g < 120.00 the transition from 
low- to high-pressure structures is characterized by an abrupt 
drop in the bond saturation, increasing from nearly full bond- 
ing (i.e., U* ~ —2.00) to a value of around —1.00. In this 
range of g, the location of the patches does not allow for the 
formation of strong bonds, thus pressure rather easily wins 
over bond saturation. 

(iii) In contrast, for g > 120.00, U* increases more 
smoothly with P* and the identified ordered structures are the 
result of a delicate tradeoff between saturation and packing. 
In particular we emphasize that for selected g values (i.e., 
g ~ 123.75 and g ~ 150.00) the respective patch decorations 
support both a high degree of bond saturation as well as the 
formation of high density lattices. As a consequence, the iden- 
tified structures are able to persist even up to high pressure 
values (i.e., P* ~ 10.00) while maintaining a relatively low 
binding energy (i.e., U* ~ —1.65). 



§ Note that for many values of the geometry parameter, the patchy interactions 
lead to - sometimes strong - deviations from perfectly symmetrical lattices. 
Therefore we refer to the structures we identified as bcc-like, fcc-like or hep- 
like to indicate such aberrations. 



The complex topology of the energy and packing fraction 
landscapes shown in Figure [2] is elucidated in the following 
by a more detailed discussion of the identified ordered equi- 
librium structures for four representative g values. Visual rep- 
resentations of the equilibrium structures at the state points 
specified in the left panel of Figure[2]by labels "a" to "h" can 
be found in the supplementary material 1 '. 

For g = 93.75 a layered structure with full bond saturation 
(£/* ~ —2.00) is formed at very low pressure (P* < 0.05): each 
layer consists of fully bonded particles, forming a honeycomb 
lattice; the inter-layer bonding is realized via the patches lo- 
cated at the north poles ("a"). Upon compression, the system 
forms a hep-like structure (i.e., AB stacking of hexagonally 
arranged layers) which is characterized by rather weak intra- 
layer (U* ~ —0.30) and stronger inter-layer (U* ~ -0.40) 
bonds ("b"). 

For the case of a regular tetrahedral patch arrangement 
(g ~ 109.47), particles arrange - as expected - at low P* in 
a bcc-like lattice ("c"). As discussed in Ref.— , the bcc-like 
configuration can be seen as two interpenetrating diamond lat- 
tices, which only interact via mutual repulsion at short dis- 
tances. The optimized lattice at T = shows another property: 
The two sublattices are slightly shifted against each other, i.e., 
particles of sublattice A do not lie exactly in the center of the 
voids of sublattice B and vice versa 0. By this rearrangement, 
optimal bond lengths can be retained at higher densities. At 
pressure values 3.20 < P* < 3.40, a variant of the bcc-like 
structure, where one of the bonds is broken in order to in- 
crease the packing fraction by nine percent, is stable. For 
larger P*, the transition into a compact fcc-like structure (with 
U* ~ -1.10 and r\ ~ 0.71, "d") takes place. We note that 
hep-like structures are also identified by the search algorithm 
as low-lying local minima at high pressure, but turn out to be 



fin Figure 2, central left panel, of the supplementary material^, the yellow 
particles would lie exactly on top of the red particles in a perfect bec lattice. 



thermodynamically unstable in this model. For a more de- 
tailed study, we refer to Ref . 

For g = 123.75, we encounter a lattice with the same bond- 
ing pattern as the aforementioned bcc/double diamond lattice 
for low pressure values ("e"). The ordered medium pressure 
phase is different from the one encountered for the regular 
tetrahedral patch decoration: particles self-assemble now in a 
rather closely packed bc-like lattice (r| ~ 0.68, "f "). However, 
the particular patch decoration allows now for a considerably 
enhanced bond saturation: bonds are formed between second- 
nearest neighbours in the bcc-like lattice, i.e., along the edges 
of the bcc-like unit cell, leading to U* ~ —1.65. This is the 
only configuration with r\ < 0.70 that is able to survive up to 
the highest pressure values investigated. 

Finally, for g = 150.00, the ordered low P* phase can be 
viewed as a stacking of staggered, hexagonally ordered double 
layers: each double layer is formed by two congruent hexago- 
nal particle arrangements, bonded to each other via the patch 
located at the north pole. The double layers themselves are 
connected via the three patches located at the basis; the bonds 
are fully saturated ("g"). For the high pressure phase (i.e., 
P* > 3.00), the system forms an ABAB... stacking of hexago- 
nally close-packed layers. Like in the low pressure configura- 
tion, the particles are essentially unbonded within these layers, 
but establish relatively strong bonds to the neighbouring lay- 
ers (with U* — —1.65, "h"). In contrast to the low pressure 
structures, the particle axes are tilted against the normal vec- 
tor of the layers, thereby breaking one bond per particle and 
considerably decreasing the interlay er distance. 

With the candidate structures suggested by the optimization 
procedure we can now proceed to the evaluation of the phase 
diagram at finite temperature (i.e. T* up to 0.3). In addition 
to the energetic and packing contributions to the thermody- 
namic potential, entropic effects come into play and config- 
urations that correspond to local enthalpy minima at T = 
might be stabilized at finite T. Thus we also consider lat- 
tices that are identified by the optimization procedure as min- 
ima that differ in enthalpy at most by twenty percent from 
the global minimum structure at some pressure value. A thor- 
ough investigation of the competition between structures cor- 
responding to local minima on the enthalpy landscape for a 
specific geometry can be found in Ref. 42 . Here, we selected 
two patch decorations, which we consider to be of particular 
interest. The first one is specified by g ~ 109.47, correspond- 
ing to the regular tetrahedral patch arrangement. As the sec- 
ond one (g = 127.50), we chose a geometry for which three 
different stable configurations have been identified at T* = 0; 
among those is the aforementioned strongly bonded be struc- 
ture, which can persist up to high pressure values. 

When the particles are decorated via a regular tetrahedral 
patch arrangement, the phase diagram shows (cf . Figure [3 
left panel) that two of the candidate structures predicted to be 
stable at T* = are able to survive over a relatively large tem- 
perature range: the bcc-like/double diamond configuration at 
low P* ("A" in FigureO and the fcc-like configuration at high 
P* ("B"). The broken bcc-like/double diamond structure on 
the other hand plays only a minor role in the temperature range 
considered: our calculations suggest that it is only stable for 
T* < 0.006. On increasing T*, an additional fcc-like phase 



("C"), which is closely related to a local minimum on the 
(T* = 0) enthalpy landscape!!, becomes stable^. Its region of 
stability emerges out of a triple point (T* = 0.064, P* = 4.39) 
and rapidly extends with increasing temperature at the cost 
of the two low temperature phases. With further increasing 
temperature, the spatially and orientationally ordered lattices 
transform via first order transitions at low pressure into the 
liquid phase and at intermediate and high pressure into an fee 
structure with perfect symmetry, where the particles are orien- 
tationally disordered (plastic crystal). At even higher temper- 
atures (not investigated here) the coexistence pressure of the 
plastic crystal and the fluid phase will gradually increase. 

Diamond cubic and diamond hexagonal structures, which 
represent local minima at very low P* values, are never stable, 
but might be stabilized by shifting the minimum of the inter- 
action potential to smaller inter-particle distances (cf. Ref.—). 
Comparing the present results with the ones from Ref. 13 , we 
note that the high density fcc-like configuration ("B"), which 
has not been taken into account for the previous calculations, 
makes the phase diagram considerably more complex. 

The phase diagram of the second system investigated (cf. 
Figure [3j right panel) displays some additional interesting 
features. With increasing T*, the predicted low- and high- 
pressure phases (bc-like, "a" and fcc-like, "Y\ respectively in 
Figure [3]) extend their region of stability at the cost of the re- 
gion of stability of the high density bcc-like phase ("p"). This 
region terminates in two very close triple points ("a"-"P"-"Y' 
and "a"-"Y'-fc c plastic crystal). At low P*, the "a" phase 
transforms into the fluid phase, while at intermediate and high 
P*, the "a" and the "Y' phases transform via first order phase 
transitions into the fee plastic crystal phase. We also observe 
a re-entrant scenario, which occurs along the coexistence line 
between the plastic crystal and the "a" phase. In a narrow tem- 
perature range (i.e., for 0.163 < T* < 0.164) the orientational 
disorder is first replaced by orientational order and is eventu- 
ally re-established upon lowering pressure. The origin of this 
behaviour is the higher compressibility of the fee plastic crys- 
tal compared to the bc-like "a" solid. The latter structure is 
stabilized by its low energy, which is only obtained when the 
bond length is close to the LJ minimum. 

The high degree of internal consistency between the two 
methodological approaches combined in this contribution is 
corroborated by the fact that the coexistence lines evaluated 
via simulations, tend towards the coexistence pressure value 
predicted by the optimization procedure for T* — >• (indicated 
by dots on the v-axes in Figure [3]) . 

The detailed structural analysis and the phase diagrams pre- 
sented in this contribution have provided quantitative evidence 
of the complex phase behaviour and the broad variety of or- 
dered equilibrium structures into which patchy particles are 
able to self-assemble. These highly complex configurations 
are the result of an intricate competition between bond for- 
mation, packing and entropy. The possibility to discriminate 
between different structural particle arrangements that are sep- 
arated only by minute energy differences in combination with 
its reliability make our conceptual approach highly suitable for 



| More specifically, the fourth-lowest local minimum predicted by the optimiza- 
tion method is structurally and orientationally very similar and immediately 
relaxes to structure "C" in the MC simulations. 




Fig. 3 (Colour online) (P*, T*) phase diagram for g ~ 109.47, corresponding to a regular tetrahedral arrangement of the patches (left) and 
g ~ 127.50 (right). Coexistence lines are calculated down to T* = 0.01. The dots on the y axes indicate the coexistence pressures at T* = 
predicted by the optimization method (neglecting the bc-broken phase); the (blue) dot in the right panel indicates re-entrancy. Labels 
correspond to the ones in Figure El right panel. 



future structural and thermodynamic investigations of patchy 
particle systems. 
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Supplementary information 



Orientational reference Hamiltonian 

When calculating free energies of systems consisting of anisotropically interacting particles via the Frenkel- 
Ladd method, the reference Hamiltonian includes an orientational term. The contribution to the entropy 
stemming from the orientational degrees of freedom can be accounted for by integrating to a reference 
Einstein crystal in which the particles are coupled to an orientational field that tends to align the particles 
with the same orientation as in the solid structure under study. It is usually convenient that the orientational 
field exhibits the same symmetry as the particles—. Therefore we used two different orientational fields for 
each of the geometries of the particles for which the phase diagram at finite T was computed. For the perfect 
tetrahedral case (with Tj symmetry), the orientation of each particle in the reference structure can specified 
by two unitary vector ao and bo that are colinear with any two specified patches, and the orientational field 
is defined as:— ^ 



where X Q is the coupling parameter, which has units of energy, and *P fljl - Qi?b,i) * s the angle formed by the 
closest patch in the instantaneous orientation of molecule i and the vector ao QV^i) i n the reference structure. 
The particles with g « 127.50 have C^ v symmetry. In this case the orientation of the particles is defined 
by a vector colinear to the patch at the "north pole" (ao) and the other is colinear to one of the other three 
patches (bo). The orientational field is the defined as: 



C/orient = £ K [sin 2 (^) + sin 2 



(1) 



i=l 




(2) 



Visual representations of selected ordered equilibrium structures at T = 

In Figures 1 to 4, we display the ordered equilibrium structures identified at T = for those systems that 
are specified in the panels of Figure 2 of the main article by labels "a" to "h". For the convenience of the 
reader, we repeat here the descriptions of these strucutres provided in the main article. 

• "a - open": layered structure with full bond saturation; each layer consists of fully bonded particles, 
forming a honeycomb lattice; the inter-layer bonding is realized via the patches located at the north 
poles 

• "b - hexagonal layers I": hep-like structure (i.e., AB stacking of hexagonally arranged layers) which is 
characterized by rather weak mfra-layer and stronger inter-layer bonds 

• "c - be I": bcc-like configuration; can be seen as two interpenetrating diamond lattices, which only 
interact via mutual repulsion at short distances; the two sublattices are slightly shifted against each 
other, i.e., particles of sublattice A do not lie exactly in the center of the voids of sublattice B and vice 
versa, as can be seen in the central left panel of figure the yellow particles would lie exactly on top 
of the red particles in a perfect bec lattice 

• "d - fc II": compact, almost close-packed fcc-like structure 

• "e - be I": lattice with the same bonding pattern as structure "c" (bcc/double diamond), but slightly 
different arrangement of the colloidal particles 

• "f - be II": rather closely packed bc-like lattice; the particular patch decoration allows for a high bond 
saturation: bonds are formed between second-nearest neighbours in the bcc-like lattice, i.e., along the 
edges of the bcc-like unit cell 

• "g - hexagonal layers III": stacking of staggered, hexagonally ordered double layers: each double layer 
is formed by two congruent hexagonal particle arrangements, bonded to each other via the patch located 
at the north pole; the double layers themselves are connected via the three patches located at the basis; 
all bonds are fully saturated 

• "h -hexagonal packed III": ABAB... stacking of hexagonally close-packed layers; particles are essen- 
tially unbonded within these layers, but establish relatively strong bonds to the neighbouring layers 




Fig. 4 Perpendicular views of the low-pressure (top) and the high-pressure (bottom) ordered equilibrium 
structure obtained for g = 93.75 (labels "a - open" and "b - hexagonal layers I" in Figure 2 of the main 
article). The colour code for the blue and green patches has been specified in Figure 1 of the main article, 
the colours red and yellow for the patchy particles have been introduced for convenience*. 



* (i) For layered honeycomb lattices and structures consisting of hexagonal layers, particles located in different layers are coloured 
in an alternating pattern, (ii) For bcc-like lattices, particles located on the vertices of the cube appear in red, while the particles 
at the center of the cube are colour yellow, (iii) In the double diamond picture, particles belonging to different non-interacting 
diamond sublattices appear in different colours, (iv) For fcc-like structures, the particles located at the vertices of the cube and 
the particles at the centers of the faces of the cube are coloured red and yellow, respectively. 
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Fig. 5 Perpendicular views of the low-pressure ordered equilibrium structure in the body-centered picture 
(top) and the double-diamond picture (center) as well as the high-pressure structure (bottom) obtained for 
g ~ 109.47 (labels "c - be I" and "d - fc II" in Figure 2 of the main article), corresponding to a regular 
tetrahedral arrangement of the patches. The colour code for the blue and green patches has been specified 
in Figure 1 of the main article, the colours red and yellow for the patchy particles have been introduced for 
convenience*. 




Fig. 6 Perpendicular views of the low-pressure (top) and the high-pressure (bottom) ordered equilibrium 
structure obtained for g = 123.75 (labels "e - be I" and "f - be II" in Figure 2 of the main article). The 
colour code for the blue and green patches has been specified in Figure 1 of the main article, the colours 
red and yellow for the patchy particles have been introduced for convenience*. 




Fig. 7 Perpendicular views of the low-pressure (top) and high-pressure (bottom), ordered equilibrium 
structure obtained for g = 150.00 (labels "g - hexagonal layers III" and "h - hexagonal packed III" in 
Figure 2 of the main article). The colour code for the blue and green patches has been specified in Figure 1 
of the main article, the colours red and yellow for the patchy particles have been introduced for 
convenience*. 



